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Abstract 

Compactly  supported  spline- wavelets  are  used 
to  process  simulated  transient  signals  in  a  noisy 
environment.  It  is  found  that  linear-phase  spline- 
wavelets  can  be  effectively  used  as  near-optimal 
time-domain  windows  to  localize  the  signal  on¬ 
sets.  This  technique  does  not  require  the  a  priori 
knowledge  of  the  signal.  Results  are  demonstrated 
numerically  by  the  detection  of  signal  arrivals  in 
noise. 


1.  Introduction 

The  arrival-time  and  spectrum  of  a  transient 
signal  contain  essential  information  about  the  po¬ 
sition  and  velocity  of  a  moving  target  in  the  Pulse 
Dopplar  mode  of  radar  operation.  Conventional 
techniques  determine  the  arrival-times  of  the  echoes 
by  either  applying  a  threshold  algorithm  or  using 
statistical  methods  after  noise  is  removed  from  the 
signal.  In  this  paper,  we  apply  the  decomposition 
algorithm  of  the  wavelet  technique  to  determine 
the  time  of  arrival  of  a  class  of  nonstationary  sig¬ 
nals  of  which  the  arrival  times  and  frequency  spec¬ 
tra  are  both  unknown. 

TVaditional  Fourier  transform  method  fails  to 
determine  the  arrival-time  since  it  produces  only 
spectral-domain  information.  Time-domain  tech¬ 
niques  without  time  and  frequency  localization  prop¬ 
erties  are  ineffective  for  this  type  of  problems.  The 
short-time  Fourier  Transform  (STFT)  generates 
time  localization  in  the  Fourier  method  by  insert¬ 
ing  a  window  function  in  the  Fourier  transform. 


The  choice  of  the  window  function  is  critical  to 
its  success.  For  high  frequency  transient  signals, 
the  time-window  must  be  very  narrow  to  produce 
good  time-domain  localization.  On  the  other  hand, 
a  wide  window  is  needed  to  deduce  the  spectral 
components  of  a  low-frequency  transient.  Because 
STFT  uses  a  fixed  window  for  all  frequencies,  it 
is  not  effective  for  wideband  signals.  The  integral 
wavelet  transform,  with  its  zoom-in  and  zoom-out 
capability,  can  be  used  in  a  wide  variety  of  sig¬ 
nal  processing  applications  [1].  In  this  paper,  we 
apply  the  compactly  supported  spline- wave  let  de¬ 
veloped  by  Chui  and  Wang  [2]  for  the  detection  of 
arrival-times  of  transient  signals. 

2.  Wavelet  series  and  integral  wavelet  transform 

Analogous  to  the  Fourier  analysis,  it  has  been 
shown  that  any  finite  energy  signal  f(t)  can  be 
represented  by  a  wavelet  series  [3] 

(1)  /(<)  = 

where 

(2)  tM0  =  2*/2tH2‘t-j) 

is  the  basic  wavelet  rp  dilated  by  a  factor  2*  and 
translated  by  j  on  the  real  line.  The  coefficients 
are  given  by  the  integral  wavelet  transform  (IWT) 
defined  by 

(3)  (W;/)(/t°>=i^/_”/<0*(4V. 

where  tp  is  the  dual  of  ip  in  the  sense  that 

(’ Pj.k.'Pt.m )  =  tj.thk.m- 


From  (3),  IWT  can  be  interpreted  as  the  pro¬ 
jection  of  the  signal  f(t)  onto  the  wavelet  space 
spanned  by  the  translates  of  the  basic  wavelet  ^ 
scaled  by  a  factor  of  a.  When  a  is  set  to  be  2~k 
and  0,j/2k,  the  wavelet  series  in  (1)  may  be  used 
to  compute  the  coefficients 

(4)  4.*  =  (W;/)  (£.£)■ 

One  advantage  in  using  wavelet  analysis  to  detect 
the  arrival-times  of  signals  is  the  variable  win¬ 
dow  size_of  offered  by  the  wavelet.  The  basic 
wavelet  ip  in  (3)  may  be  regarded  as  a  window 
function  whose  width  is  flexibly  controlled  by  the 
scale  value  a  to  produce  the  zoom-in  and  zoom- 
out  effects  used  for  examining  a  signal.  If  a  is 
chosen  to  be  large,  the  time-window  is  wide  and 
the  IWT  produces  good  low-frequency  resolution. 
Conversely,  a  small  value  of  a  generates  a  nar¬ 
row  iime-window  for  time-domain  indication  of 
the  arrival-time  of  a  signal. 

3.  Spline-wavelets  and  algorithms 

The  compactly  supported  spline- wave-lets  de¬ 
veloped  by  Chui  and  Wang  [2]  are  explicitly  rep¬ 
resented  by  a  linear  combination  of  the  5-spline 
functions.  If  Nm{t)  denotes  the  cardinal  5-spline 
of  order  n,  then  the  mth  order  compact  supported 
spline  wavelet  is  given  by 

2  m  — 2 

<M*)  =  ^rr  E  (-n^mO'  +  D 

J= 0 

(5) 

The  cubic  (m  =  4)  spline- wavelet  and  its  Fourier 
transform  are  illustrated  in  Figures  1  and  2,  while 
the  signal  processing  properties  of  the  spline-wavelets 
for  m  =  2  and  4  have  been  evaluated  and  tabu¬ 
lated  in  Tables  I  and  2.  It  can  be  seen  from  Table  1 
that  the  area  of  the  rectangle  (Arpt  x  Aipw)  on  the 
time-frequency  plane  is  approximately  equal  to  0.5 
for  the  cubic  spline- wavelet  which  is  the  lower  limit 
due  to  the  uncertainty  principle.  Numerically,  we 
can  show  that  as  the  order  of  the  spline  wavelet 
is  increased,  the  area  of  the  rectangle  in  the  time- 
frequency  plane  approaches  the  lower  limit.  Hence 
the  window  property  of  the  spline- wavelets  is  near 
optimal. 

It  is  also  easy  to  see  from  Figure  2  that  the 
spline-wavelet  is  the  impulse  response  of  a  band¬ 
pass  filter  The  IWT  can  be  interpreted  as  a  band¬ 
pass  filtering  of  the  signal  f{t)  where  the  center 


frequency  and  the  bandwidth  is  controlled  by  the 
parameter  a.  For  the  cubic  spline- wavelet,  the 
first  sidelobe  of  the  bandpass  filter  is  better  than 
-90  db. 

Furthermore,  the  geometric  symmetry  of  the 
spline- wavelets  guarantees  the  linear  phase  prop¬ 
erty.  Consequently,  the  spline-wavelets  may  be  re¬ 
garded  as  linear  phase  bandpass  FIE  filters.  Two 
algorithms  have  been  constructed  using  spline  wavelets 
with  these  properties  for  processing  transient  sig¬ 
nals  in  real-time. 

The  first  is  a  read-time  algorithm  for  project¬ 
ing  the  sampled  signal  f(n)  onto  the  spline  space. 

The  goal  here  is  to  find  the  coefficients  of  a  spline 
series  approximating  the  signal.  The  projection 
formula 


(6)  m  *  /.(<) 

=  E  {e&U-A*)} Nm  (h -  J) 

with  h0  =  tmh,  has  been  used  in  [4].  For  cubic 
spline  (m  =  4),  the  weight  coefficients  are 
tabulated  in  [4]. 

The  second  algorithm  is  to  decompose  the  ap¬ 
proximating  signad  fa(t)  into  several  orthogonal 
components  at  various  resolutions.  Assuming  that 
fa(t)  is  represented  at  the  tth  resolution,  (6)  can 
be  rewritten  into  a  spline  series 

(7)  /«(<)  =  Ec‘^m(j  -;)=/*(<) 

which  can  be  decomposed  into  wavelet  series  at 
lower  resolutions 

(8)  fk  =  9k-l  ©  9k-2  ©  ’  •  •  ©0*-l  +  fk-L 

where  L  is  an  arbitrary  positive  integer  and  ® 
denotes  orthogonal  sums.  The  wavelet  series  for 

MO  is 

(9)  MO  =  E  (^2 F/i  “  • 

where  the  wavelet  coefficients  dp j  are  obtained  by 
a  moving  average  operation 

(10)  dp.j  =  E  bn-2jcp+l,n 


'  i '■  C-c.t-3 


□  □ 
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followed  by  subsampling  and  the  spline  coefficients 
cpj  are  also  obtained  by  a  similar  operation,  namely: 

(11)  Cp  j  —  "y  ^  an-2jCp+].,n- 

n 

The  coefficient  sequences  {at}  and  {6t}  are  tabu¬ 
lated  in  [1]. 

4.  Numerical  examples 

A  wide  variety  of  physical  phenomena  can  be 
represented  by  sinusoidal  signals  bounded  by  ex¬ 
ponential  decay  enve¬ 

lopes.  A  signal  consisting  of  three  different  fre¬ 
quencies  is  given  by: 

3 

(12)  y(t)  =XVA,(",,VW,‘<'(<-fO 

i=l 

where  A*  >  0  are  decay  constants  of  the  envelopes; 
f,  are  the  arrival-times;  and  w,-  are  the  signal  fre¬ 
quencies.  This  type  of  signals  has  been  studied 
in  [5]  using  a  one-sided  exponential  window  which 
matches  the  signal  waveforms.  We  analyze  these 
signals  by  the  cubic  spline-wavelet  using  the  two 
algorithms  outlined  in  Section  3.  In  our  exam¬ 
ple,  the  signal  is  corrupted  with  Gaussian  noise. 
The  times  of  arrival  of  the  three  exponentials  are 
clearly  shown  in  Figures  3b  and  3c. 

We  also  use  two  overlapping  modulated  sinu¬ 
soids  of  the  same  frequency  to  simulate  a  receiver 
pickup  of  the  direct  signal  and  its  echo.  The  echo 
amplitude  is  5%  of  the  direct  signal.  Figure  4a 
displays  the  received  signal  (direct  +  echo)  and 
the  echo  alone.  Figure  4b  shows  the  echo  sepa¬ 
rated  from  the  received  signal  with  its  arrival  time 
clearly  indicated. 

5.  Conclusions 

We  have  demonstrated  a  technique  to  detect 
the  arrival-times  of  transient  signals  using  com¬ 
pactly  supported  spline- wavelets.  Unlike  other  tech¬ 
niques  which  depend  on  match  filtering  or  statisti¬ 
cal  algorithms,  our  algorithms  decompose  the  sig¬ 
nal  into  various  frequency  bands  to  apply  thresh¬ 
olding  techniques.  The  results  of  our  numerical 
experiments  agree  very  well  with  the  designed  sig¬ 
nals. 
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Figure  1  Cubic  spline-wavelet  rp 


Figure  2  Fourier  transform  of  rp 


Figure  3a  Original  signal  with  noise 
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Figure  3c  Onset  of  the  low-frequency  signal 


Wavelet  Gas 


B- wave  let 


Dual- wavelet 


m=2 

UA\ 

0.390 

5.406 

2.493 

m=4 

i im 

0.542 

5.164 

0.932 

m=2 

TO 

4.778 

2.350 

iiiHim 

0.S73 

4.172 

1.041 

Table  l :  Centers  and  widths  of  spline- 
windows 


Figure  3b  Onsets  of  the  inid-  and  high 
frequency  signals 


Figure  4 


